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ABSTRACT 

Sliding friction between crystal lattices and the physics 
of cold ion traps are so far non-overlapping fields. Two 
sliding lattices may either stick and show static friction 
or slip with dynamic friction; cold ions are known to form 
static chains, helices, or clusters, depending on trapping 
conditions. Here we show, based on simulations, that 
much could be learnt about friction by sliding, via e.g. 
an electric field, the trapped ion chains over a periodic 
corrugated potential. Unlike infinite chains where, ac- 
cording to theory, the classic Aubry transition to free 
sliding may take place, trapped chains are always pinned. 
Nonetheless we find that a properly defined static fric- 
tion still vanishes Aubry-like at a symmetric-asymmetric 
structural transition, ubiquitous for decreasing corruga- 
tion in both straight and zig-zag trapped chains. Dy- 
namic friction can also be addressed by ringdown oscilla- 
tions of the ion trap. Long theorized static and dynamic 
one dimensional friction phenomena could thus become 
exquisitely accessible in future cold ion tribology. 




FIG. 1: (a) and (b) sketch of proposed ion trap configuration 
in a periodic potential; (c) sketch of the straight and zig-zag 
chain equilibrium geometries determined by different aspect 
ratios of the trapping quadrupolar potential. The lateral zig- 
zag amplitude has been artificially magnified. 



INTRODUCTION 

The field of sliding friction has been recently revived 
thanks to experimental and theoretical advances espe- 
cially connected with nanosystems, with brand new op- 
portunities to grasp the underlying phenomena of this 
important and technologically relevant area. Unlike 
macroscopic friction, where contact irregularities domi- 
nate, the nanoscale offers perfectly defined crystal lattice 
facets that can mutually slide. One fundamental piece of 
physics arising in this limit is the possibility of frictionless 
sliding when two perfect lattices are mutually incommen- 
surate - their periodicities unrelated through any ratio- 
nal number- an idealized situation sometimes referred to 
as superlubric [1 j. Graphite flakes sliding on graphite 
were for example shown to be pinned by static friction 
- the finite threshold force F s necessary to provoke mo- 
tion - when aligned, but to turn essentially free sliding 
once a rotation makes them incommensurate [2 . An 
additional prediction, rigorously proven for one dimen- 
sional (ID) infinite harmonic chains sliding in a periodic 
potential (the "corrugation" potential, in frictional lan- 
guage), is that the incommensurate vanishing of static 
friction will only occur so long as the sliding chain is 
sufficiently hard compared to the corrugation strength. 



When instead that strength exceeds a critical value, the 
chain becomes locked, or "pinned", to the corrugation, 
thus developing static friction despite incommensurabil- 
ity. Unlike most conventional phase transitions, this cele- 
brated Aubry transition [3j[4] involves no structural order 
parameter and no breaking of symmetry - just a change 
of phase space between the two states, the unpinned and 
the pinned one. The concept of pinned and unpinned 
states is by now qualitatively established in the sliding of 
real incommensurate 3D crystal surfaces [5 . Yet, the on- 
set of static friction in the prototypical ID chain, where 
the Aubry transition is mathematically established, has 
never been experimentally validated. Even less is known 
about dynamic friction, about which we have no other 
insight than generic linear response [6], and some data 
on 3D rare gas overlayers [7 both suggesting a quadratic 
increase with corrugation. 

Cold atomic ions can form linear chains. Despite their 
mutual Coulomb repulsion, ions can be corralled inside 
effective potential traps generated by quadrupolar radio- 
frequency electrodes [8]. The low temperature ion equi- 
librium geometry is determined by the aspect ratio asym- 
metry of the confining potential, according to the effec- 
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tive hamiltonian [9j [10] 
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where = (x^ y%,Zi) is the position of the z-th ion, m its 
mass, e its charge, and the three terms represent kinetic 
energy, Coulomb repulsion, and the effective parabolic 
trapping potential respectively, with transverse and 
longitudinal trapping frequencies u± and cjy. As is 
known experimentally [IT] and theoretically p~2j [13] , by 
decreasing the asymmetry aspect ratio R = oj\\/oj± the 
trap potential deforms from spherical to cigar-shaped, 
and the equilibrium ion geometry changes from 3D 
clusters to helices to chains, initially zig-zag and finally 
linear and straight along z, as in Fig. [ljc). As many 
as a hundred ions may be stabilized in a linear con- 
figuration, with a few jam ion- ion spacing ao, typical 
repulsion energy e 2 /ao ~ 0.5 K, and temperatures below 
~ lfiK. So far, cold ion chains raised interest in view 
of promising applications for spectroscopy [T4tfl6] . fre- 
quency standards [T7J [18] , study and control of chemical 
reactions [19] , and quantum information [20-23 . We are 
concerned here with the possibility that they could be 
of use in the field of friction. An interesting observation 
was made by Garcia-Mata et al. [24 , that ion chains 
in an additional incommensurate periodic potential 
(produced, e.g., by suitable laser beams) would, if only 
the spatial inhomogeneity of the trapping potential 
could be neglected, resemble precisely the ID system 
where an Aubry transition is expected as a function 
of corrugation strength. Unfortunately, the trapping 
potential itself introduces an additional and brutal break 
of translational invar iance, at first sight spoiling this 
neat idea. The ion-ion spacing is only uniform near the 
chain center, increasing strongly towards the two chain 
ends, where it diverges owing to the trap confinement 
(see e.g., Fig. [I]). Incommensurability of the chain with 
the periodic potential is therefore lost, and seemingly 
with that all possibilities to study phenomena such as 
the vanishing of static friction at the Aubry transition, 
a transition known to depend critically on a precisely 
defined value of incommensurability [25] . 

However, this conclusion is overly pessimistic. In con- 
texts unrelated to friction, a finite chain in a periodic 
potential [26H28] is known to turn the would-be Aubry 
transition to a structural phase transition. The transition 
survives, as we will show, despite the addition of the con- 
fining trap potential. The distorted state, prevailing at 
strong corrugation, exhibits under suitable conditions a 
change of symmetry and an accompanying displacive or- 
der parameter, measuring the breaking of inversion sym- 
metry about the chain center. Simulations carried out by 
pulling the trapped chain with a longitudinal force (such 
as in principle provided by a static external electric field) 
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FIG. 2: Zero temperature ion chain equilibrium center-of- 
mass coordinate 5zcm versus periodic corrugation amplitude 
Uo for the linear chain (red) and zig-zag chain (blue). Posi- 
tion, time and energies in dimensionless units (see methods). 
The insets zoom on the linear chain and the zig-zag central 
portions, highlighting the symmetry breaking taking place at 
Uq = U c and above. 



reveal that the distortion transition is associated with a 
continuous onset of static friction from zero to finite, thus 
mimicking at finite size, and despite the trapping poten- 
tial, the static friction onset theoretically expected at the 
ideal Aubry transition. 

In the rest of this work explicitly simulated gedanken 
experiments will demonstrate that trapped cold ions may 
in fact represent an idesal system for future nanotribo- 
logical studies addressing both static and dynamic fric- 
tion. In static friction, the celebrated Aubry transition 
will become experimentally accessible for the first time. 
In dynamics, the friction rise with increasing corrugation 
should become measurable giving new impulse to future 
theory work. 



RESULTS 

Model and symmetry breaking transition. 

We simulate by means of standard (weakly damped) 
classical molecular dynamics N = 101 ions trapped by 
the parabolic potential of Eq.Q, and seek first the static 
(T = 0) equilibrium geometry. Below a critical aspect ra- 
tio Ro = 1.18 • 10 -3 the equilibrium geometry is a linear 
chain, while at R = Ro there is a transition to a pla- 
nar partly zig-zag structure (see Fig. [I] and Supplemen- 
tary Movie 1), further evolving to a helical structure, and 
eventually to a 3D cluster when R 1 [12]. Choosing 
R = 0.5 • 10 -3 for the frictional simulation of ID trapped 
linear chain, we introduce an additional periodic "cor- 
rugation" potential U = Z7 cos(^ L z). To address the 
best studied case of Aubry transition, we choose the cor- 
rugation wavelength A to be "golden ratio" incommen- 



3 



surate relative to the chain center ion- ion spacing ao, 
A = 2ao /(l + \/5) [3j [4] . Preserving the chain's inversion 
symmetry, the periodic potential phase is chosen so that 
the central ion sits at a potential maximum. At van- 
ishing external force, the chain equilibrium configuration 
can then be followed while adiabatically increasing the 
corrugation amplitude Uo- As shown in Fig. |2j the chain 
ground state remains fully left-right symmetric right up 
to a critical corrugation amplitude U c = 0.01627, (in di- 
mensionless units, see Methods), where it sharply turns 
asymmetric, the center of mass (CM) moving right or 
left by an amount Szcm(Uo) (see Supplementary Movie 
2). This transition, similar but not identical to that 
known for free finite chains [26] [28] is the remnant of 
the Aubry transition of the ideal infinite chain. In both 
the finite and the infinite chain, the probability to find an 
ion (here the central ion) sitting exactly at a periodic po- 
tential maximum vanishes above the critical corrugation 
strength, where the ion moves out to a minimum distance 
ip away from the potential maximum [4, 29 . Unlike the 
Aubry transition that has no structural order parameter, 
the (T = 0) finite chain transition is structural, char- 
acterized by a displacive order parameter Szcm, Ising- 
like since the central ion can identically fall left or right 
of the centre. This order parameter grows continuously 
as \Uo — U C \P at the transition, (see Fig. [2J. We find 
j3 — 0.496 ±0.005 for the trapped linear ion chain - prac- 
tically coincident with the exponent k = 0.5 known for 
finite chains with short range forces and without trap- 
ping potential [28] . Moreover, upon reducing the confin- 
ing shape asymmetry by increasing the aspect ratio to 
R = 1.2 • 10 -3 , the chain equilibrium geometry evolves 
from linear to a partly planar zig-zag (Fig. [2| , still sym- 
metric about z = Oat[/o = 0. For increasing corru- 
gation a symmetry-breaking transition similar to that of 
the linear chain takes place here too at U c = 0.01357 - 
a smaller value since in the zig-zag the ions are further 
apart - and for such case /3 = 0.504 ± 0.005 (see Supple- 
mentary Movie 3). For both the linear and the zig-zag 
chain the critical amplitude U c also depends upon the 
incommensurability ratio A/ao (not shown). While these 
results are valid at T=0, it should be stressed that as in 
all small size systems, temperature will cancel all sym- 
metry breaking and order parameters, owing to thermal 
jumps taking place above the energy barrier A which sep- 
arates the two opposite Ising like order parameter valleys. 
However, the mean time lapse between jumps diverges 
exponentially when T « T5 ~ A/&b, the trap blocking 
temperature. Except too close to U c , ion traps are rou- 
tinely cooled down to temperatures of order 10 -7 (in our 
units e 2 /d ~ 6.3 K) way below the blocking temperature 
(for example T5 ~ 180/iK even for Uo = 0.014, barely 
above the critical U c = 0.01357). Even at finite temper- 
atures therefore, the ion chain will still break inversion 
symmetry and exhibit a distortive order parameter for 
macroscopically long times similar to the T = results 



detailed above. 



Determination of the static friction force. 

Thus far the description of the ion trap in a periodic 
corrugation, while still free of external forces. To study 
static friction we apply a uniform force F, such as that 
of an electric field parallel to the chain axis. The usual 
procedure used to determine static friction - the small- 
est force capable of causing a pinned-to-sliding transition, 
following which the chain is depinned with unlimited slid- 
ing - does not work here, since the field-induced cen- 
ter of mass displacement Szcm(F) is finite by necessity, 
owing to the confining trap potential. We circumvent 
that problem by monitoring another quantity, namely 
the restoring force Fr that must be applied to the chain 
in order to shift its CM coordinate back exactly to zero 
- keeping the chain's central ion exactly on top of the 
corrugation maximum at z = 0. In the infinite free chain 
this alternative definition of static friction coincides with 
the standard pinned-to-sliding Peierls-Nabarro threshold 
force Fpn [4j [30] , identifiable with the static friction F s 
of the system. In the pinned Aubry state the critical force 
Fp n is exactly such as to drive to zero the minimum equi- 
librium deviation 1/; of one "kink" away from the near- 
est corrugation maximum, eventually triggering the free 
motion which initiates the pinned-to-sliding transition of 
the whole system [4j [31] . Analogously, for trapped ions 
a growing external force F gradually causes the central 
atom deviation from the potential maximum to decrease, 
and eventually to vanish at F = Fr. We conclude that 
Fr measures the static friction even in a trap, where 
proper free sliding is impeded. As shown in Fig j3] (inset), 
for Uo < U c the central ion of the trap spontaneously 
remains precisely at zero, whence Fr = 0. However at 
Uo > U c where the force- free displacive order parame- 
ter Szcm turns spontaneously nonzero, a finite restor- 
ing force Fr must be applied in order to cancel it. The 
growth of the symmetry restoring force Fr with Uo ( 
Fig. [3| is, despite the larger threshold, generally close to 
the Aubry static friction expected theoretically for the in- 
finite ideal chain [IJ [3] , when calculated for same chain 
parameters ao, Uo and golden ratio A, as those at the 
trapped chain center. We have thus established a con- 
nection between the trapped chain restoring force Fr and 
the Peierls-Nabarro static friction of the infinite chain, 
of tribological significance pQ. Numerical differences be- 
tween the two are of course hardly surprising. The finite 
size Nao cuts off the Aubry critical region where the cor- 
relation length exceeds this size. For the same reason the 
chain depinning exponent x, Fr ~ \Uo — U c \ x , which we 
find to be about 1.85 ± 0.2 differs from the higher ex- 
ponent ~ 3 expected for the infinite free chain for short 
range forces |4j [29] [30] . The effective static friction of 
linear and zig-zag chains also differ slightly (see Fig. [3| , 





FIG. 3: Effective static friction measured by the restoring 
force Fr necessary to symmetrize the system about the cen- 
ter, versus corrugation amplitude Uo for the trapped linear 
ion chain (red) and zig-zag chain (blue). The insets illustrate 
the procedure used. Black dashed line: calculated theoreti- 
cal static friction at the Aubry transition of the ideal infinite 
golden ratio incommensurate chain. Arrows mark the vanish- 
ing point of effective static friction in each case. 
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FIG. 4: Oscillatory kinetic ringdowns for the linear (upper 
panel) and the zig-zag (lower panel) chains. The two panels 
show the zcm coordinate as a function of time for different 
values of Uo < U c \ red Uo = 0, green Uo = 0.0001, blue 
Uo = 0.0002, magenta U = 0.0003 and cyan Uo = 0.0004. 



the softer zig-zag chain easier to pin, owing to its larger 
compliance than the linear chain. Overall, both cases 
provide a measurable realization of the sharp onset of 
static friction at the Aubry transition, valid beyond the 
limitations of a trapped chain. 



Dynamic friction. 

Besides static friction and its Aubry-like vanishing 
transition, ion traps could in addition provide nontriv- 
ial insight on dynamic friction. The onset of dynamical 
friction in depinned infinite chains is a very subtle event 
as discussed in literature [32] ; the pinned-to-sliding onset 
has been simulated and elegantly described by Braun and 
coworkers [4j [31] . In tribology we are primarily interested 
in the growth, for incommensurate sliders at finite sliding 
speed, of dynamical friction with increasing corrugation. 

To simulate the chain dynamical behaviour we initially 
displace the ions away from their equilibrium configura- 
tion with a static field, then remove the field and follow 
the spontaneous "ringdown" damping of the bodily os- 
cillation of the whole chain in the trap potential. For an 
initial large CM displacement of order 10 ao the chain os- 
cillation damping is readily observable, and can for weak 
corrugation Uo <C U c be measured even in the absence of 
an external heat bath. The initial potential energy (here 
about 0.1) is converted to Joule heating of the chain, at a 
ringdown rate which is a direct measure of dynamic fric- 
tion. Actually, a trapped chain has an intrinsic dynamic 
friction even in the absence of corrugation, because of the 
trap-potential-induced transfer of kinetic energy from the 
CM motion to internal chain phonon modes. However, 



this effect is small; the conversion and the corresponding 
dynamic friction is strongly enhanced by the corruga- 
tion potential, as shown in Fig. [4] and in Supplementary 
Movies 4 and 5 (for the linear and zig-zag configurations 
of the chain, respectively). Simulation of the chain oscil- 
lations yields a ringdown amplitude decay for both linear 
and zig-zag chains, of the form Aoexp — [t/r]. We find 
that the lifetime t(Uq) strongly decreases for increasing 
Uo. The corresponding increase of dynamic friction Fr>, 
directly proportional to r _1 , is to a good approximation 
quadratic, F& ~ a + bU§ for weak corrugation (Fig. [5|. 
A quadratic increase is just what the linear response the- 
ory predicts based on Fermi's golden rule [6j[33]. It also 
agrees with experimental Quartz Crystal Microbalance 
monolayer friction data of Coffey and Krim [7]. There 
is in the literature no dynamic friction theory for the 
infinite linear chain to be used for comparison with the 
trapped chain results of Fig. [5j the present results can 
therefore be considered a first exploratory result in this 
direction. The exceedingly small values of a speak for 
the strong harmonicity of the chain, whose CM motion 
can only with great difficulty excite the higher frequency 
intra-chain phonons in the absence of corrugation. The 
quadratic increase with Uo reflects the effective open- 
ing of anharmonic excitation channels through Umklapp 
scattering. Larger oscillation amplitudes could not be 
readily simulated because the Joule heat liberated would 
in that case be excessive and would destroy the chain 
in the lack of some external dissipation mechanism. A 
fully nonlinear behavior of dynamic friction would be ex- 
pected to emerge in that case, with outcomes strongly 
dependent upon the heat disposal mechanism. 
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FIG. 5: (a) Inverse decay time r _1 of CM ringdown oscilla- 
tions, measuring the dynamic friction of the ion chain, linear 
(red) and zig-zag (blue). The inset illustrates the procedure 
used to obtain it. Continuos lines represent fits of the form 
a + bUo (a — 0.014 and b — 867.82 for the linear chain and 
a — 0.0 and b — 250.31 for the zig-zag). 



DISCUSSION AND CONCLUSIONS 

The remarkable static and dynamic friction properties 
emerging for trapped ion chains seem well amenable to 
experimental study. The exceedingly small values of the 
distortion 5zcm necessary to establish an accurate value 
for the effective static friction Fr will require a fine tun- 
ing of corrugation and applied force. Upon increasing 
the force however, the instant F = Fr is reached and 
the symmetric position Szcm = is attained, the chain 
will spontaneously jump over the barrier to the oppo- 
site Ising order parameter valley. By e.g., alternating in 
time rightward and leftward forces, the effective static 
friction Fr could become readily detectable in experi- 
ment as the threshold force magnitude for the chain's 
CM to sway visibly leftward and rightward. The present 
static friction study is conducted at zero temperature, 
ignoring quantum effects, without external sources of ei- 
ther heating, cooling or damping, and assuming that a 
sufficiently strong effective periodic corrugation poten- 
tial can be laser generated. While each of these points 
requires careful consideration, none, as will be discussed 
below, is in our view fatal to the results presented for 
static friction. The proposed dynamic friction experi- 
ment would heat up the trapped ion chains, and cause 
their destruction once the Joule heating is excessive. We 
verified however that even at a temperature of ~ 0.55 
K, orders of magnitude higher than the effective chain 
Debye temperature of ~ 145/iK (considering Mg + ions 
separated by ao = 5/im), the chain is still stable in the 
strong trapping potential. 

In conclusion, the field of cold atoms has had a tremen- 
dous development, recently covering important interdis- 
ciplinary areas such as Anderson localization [34]. Now 



the potential tribological impact of cold ion results could 
also be significant. The onset of static friction is a rich 
problem of current interest [35j [36] , and so is the dynamic 
friction growth with corrugation [7 . Most importantly, 
ion chains represent a new clean design, free of most com- 
plications of real crystal lattices and with many more con- 
trollable parameters for a study of basic nanofrictional 
properties. 

METHODS 

Simulation Details. The equations of motion de- 
rived from (fiT) turn dimensionless by defining the quan- 
tity d 3 = e /muj\, and expressing distance in units of 

time in units of 1/uj± and the mass in units of m. 
With this definition, energy and force are also dimension- 
less in units of e 2 /d and e 2 /d 2 respectively. Simulations 
are carried out using classical molecular dynamics, when 
necessary with a weak frictional damping, and mostly 
at zero or low temperature. The equations of motion 
have been integrated using a velo city- Ver let algorithm 
with a time step At = 5 • 10 -3 time units. The thresh- 
old trap potential asymmetry Ro for linear /zig-zag struc- 
ture transition is determined with adiabatic increments 
AR = 10~ 4 starting from R = 0.5 • 1(T 3 . The criti- 
cal periodic corrugation potential amplitude U c for the 
symmetric-asymmetric transition of linear and zig-zag 
chains is obtained by adiabatic increments AUo = 10 -5 
starting from an initial value Uo = 0.012 for the zig-zag 
configuration and Uq = 0.015 for the linear chain. In 
all the adiabatic increment procedures a small viscous 
damping — 7^, with 7 = 0.01, is applied to all ions in 
order to eliminate the excess energy. With our choice of 
7 a total time of 2500 between successive increments is 
enough to damp all the excess energy away. 

Physical parameters and orders of magnitude 

a) Temperature and quantum effects. The present 
study is conducted conventionally at zero temperature, 
ignoring quantum effects, without external sources of 
heating, cooling or damping, and assuming that a suf- 
ficiently strong periodic potential can be laser generated 
for the transition to take place at the given incommensu- 
rability provided by the laser field. Even though temper- 
ature is in reality finite (typically of order 100 nK), the 
previously reported comparison between left and right 
chain distortion energies with ksT implies a symmetry- 
broken state lifetime much longer than any measurement 
time, justifiying the zero temperature treatment, where 
data and their significance are clearer. In this regime 
moreover, quantum zero-point effects could in principle 
alter some of the details close to the depinning thresh- 
old, but are otherwise not expected to affect the overall 
friction picture, both static and dynamic. 

b) Corrugation potential magnitude. To produce the 
effective periodic corrugation potential, the wavelength 
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£ of a laser light must fit an electronic excitation in the 
spectrum of the chosen ion species (e.g. 280 nm for Mg + , 
729 nm for Ca + or 399 nm for Yb + ), with a typical 
ao = 5/im the golden mean incommensurate periodic po- 
tential wavelength is A = 3000nm ^> £, nevertheless the 
desired A value for the effective periodic potential can be 
obtained by laser beam crossing at a chosen angle [24]. 
Experimentally, it is mandatory to realize a sufficient cor- 
rugation amplitude U = Uq for the Aubry transition to 
occur. Assuming for example Mg + ions and golden ra- 
tio incommensurability, the critical amplitude of the ef- 
fective periodic potential U c requires a laser intensity of 
some KW/m 2 , which seems entirely within reach. On the 
other hand, the Aubry pinning transition becomes, as is 
well known, easier with any chain-corrugation incommen- 
surability different from the golden ratio [25] , whereby in 
general the transition would occur with even with weaker 
corrugation amplitudes than Uq. 

c) Phase misalignment effects. We assumed so far the 
possibility to position the corrugation potential maxi- 
mum precisely at the minimum of the harmonic trap, a 
feat which may be experimentally difficult. Slight shifts 
of the corrugation potential phase away from that point 
would however simply smear the Aubry pinning onset in 
the same way a weak symmetry breaking field transition 
would affect a phase transition. 

d) Chain heating during ringdown. In the oscillation 
ringdown simulations, the initial potential energy differ- 
ence 8.9, due to the chain displacement of 10 ao, is grad- 
ually converted to heating of the chain. With the mass of 
Mg + , the final temperature reached is approximately 0.55 
K, orders of magnitude higher than the effective chain 
Debye temperature of ~ 145/iK. Despite that, we ver- 
ified by direct simulation that at this temperature the 
chain is still stable in the strong trapping potential. 
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